MAGNETIC-FIELD MEASUREMENTS OF T TAURI STARS 
IN THE ORION NEBULA CLUSTER 1 



Hao Yang 2 

Department of Physics & Astronomy, Rice University, 6100 Main St. MS- 108, Houston, 

TX 77005 

haoyang@colorado . edu 

Christopher M. Johns-Krull 
Department of Physics & Astronomy, Rice University, 6100 Main St. MS-108, Houston, 

TX 77005 

cmj@rice . edu 

Received ; accepted 

DRAFT 



1 Based on observations obtained at the Gemini Observatory, which is operated by the 
Association of Universities for Research in Astronomy, Inc., under a cooperative agreement 
with the NSF on behalf of the Gemini partnership: the National Science Foundation (United 
States), the Science and Technology Facilities Council (United Kingdom), the National Re- 
search Council (Canada), CONICYT (Chile), the Australian Research Council (Australia), 
Ministerio da Ciencia e Tecnologia (Brazil) and SECYT (Argentina). 

2 Now at JILA, University of Colorado, 440 UCB, Boulder, CO 80309-0440. 



- 2 - 



ABSTRACT 



We present an analysis of high-resolution (R ~ 50, 000) infrared K-band 
echelle spectra of 14 T Tauri stars in the Orion Nebula Cluster. We model Zeeman 
broadening in three magnetically sensitive Ti I lines near 2.2 /iin and consistently 
detect kilogauss-level magnetic fields in the stellar photospheres. The data are 
consistent in each case with the entire stellar surface being covered with magnetic 
fields, suggesting that magnetic pressure likely dominates over gas pressure in the 
photospheres of these stars. These very strong magnetic fields might themselves 
be responsible for the underproduction of X-ray emission of T Tauri stars relative 
to what is expected based on main-sequence star calibrations. We combine these 
results with previous measurements of 14 stars in Taurus and 5 stars in the TW 
Hydrae association to study the potential variation of magnetic-field properties 
during the first 10 million years of stellar evolution, finding a steady decline in 
total magnetic flux with age. 

Subject headings: open clusters and associations: individual (Orion) — infrared: 
stars — stars: magnetic fields — stars: pre-main sequence 
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INTRODUCTION 



T Tauri stars (TTSs) are a class of low-mass variable stars. Typically a few million 
years (Myr) old, TTSs have recently formed out of dense molecular cloud cores and are 
evolving toward the main sequence along their Hayashi tracks. These young solar-type stars 
display many spectral pecularities, often including excess continuum emission at ultraviolet 
through infrared (IR) wavelengths, strong and variable Ha and Ca II H and K emission 
lines, both red- and blue-shifted absorption features in Balmer lines, and forbidden line 
emission. These spectral characteristics are most prominent in the spectra of so-called 
classical TTSs (CTTSs), which are TTSs still surrounded by dusty accretion disks. Another 
category of TTSs, naked TTSs (NTTSs), do not appear to have such disks around them 
and show little or no sign of accr etion. For a comprehensiv e int roduction to T TSs, readers 



are directed to recent reviews by 



Menard fc Bertout 



(119991 ) and lPetrovi (]2003h . 



Stellar magnetic fields are thought to play an important role during the evolution of 
TTSs. Strong magnetic fields are b elieved to directly co ntrol the interaction between CTTSs 



and their circumstellar disks (e.g., 



Bouvier et al 



20071 ). Magnetic fields are thought to 



truncate the disks at or near th e corotation radius, forcing disk material to accrete along the 



field l ines onto the central star (ICamenzind 



1993 



Shu et al. 



1994 



Paatz &: Camenzind 



support of this picture ( iBouvier et al. 



1990 



Konigl 



1991 



Collier Cameron fc Campbell 



19961 ) . While there is abundant evidence in 



20071 ) . observational detection of magnetic fields is 



difficult, and there remain relatively few magnetic field measurements on young stars. 

The most successful approach thus far for measuring the total magnetic flux on cool stars 



has been measuring the Zeeman broadening of spectral lines in unpolarized light ([R obinson 



Saar 


1988 




Valenti et al. 


1995: 


Johns-Kru 


11 k Valenti 


1996 


Yang et al. 


2005 


> 


2008 




Johns-Krull 


2007 


). For any given 



Johns-Krull et al. 



1999b . 
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splitting resulting from the magnetic field is 



AA [mi 



A7rm e c 2 



\ 2 gB = 4.67 x 10- { \ 2 gB , 



where g is the Lande ^-factor of the transition, B is the strength of the ma gnetic field in kG 



Johns-Krull et al. 



(11999bl hereafter Paper 



and A is the wavelength of the transition in A. 
I) detected Zeeman broadening of the Ti I line at 2.2233 fim on the CTTS BP Tau and 
obtained a field strength of B = 2.6 ± 0.3 kG av eraged over the en t ire su rface of the star. 



Johns-Krull et al. 



(120041 hereafter Paper 



From the analysis of four Ti I lines near 2.2 /im, 
II) measured a strong magnetic field of 2.5 ± 0.2 kG on the NTTS Hubble 4. In addition 
to studying the Ti I lines, Paper II also examined several CO lines near 2.3 /mi which have 
substantially reduced magnetic sensitivity and found no excess broadening beyond that due 
to stellar rotational and instrumental broadening. 

Consistent detection of strong fields on TTSs raises the questio n of the orig i n of 



Johns-Krulll (120071 ) 



these fields. In the largest study to date of magnetic fields on TTSs, 
used a sample of 14 stars to look for correlations between the measured magnetic-field 
properties and stellar propert ies that might be i mportant for dynamo action. No significant 



correlations were found, and 



Johns-Krulll (120071 ) speculated that the fields seen in these 



youn g stars 



[e.g., 



Tavler 



may 



De entrainec 



1987 



Moss 



20031 ) 



interstellar fields le: 



Yang et al 



(12005 



t oyer from the star formation process 



20081 . hereafter Papers III and IV) 



studied six TTSs in 



the T W Hydrae association (TWA). Roughly 10 Myr old (e.g. 



Barrado Y Navascues 



200 6]), the TWA stars are substantially older than the Taurus stars, 



whose age is ~ 2 Myr (jPalla fc Stahler 



20001 ) . Paper IV reported that the total unsigned 



magnetic flux of the stars appears to decrease with time from 2 Myr to 10 Myr, which 
further supports a primordial origin of the magnetic fields on TTSs. 

In this work, we investigate the magnetic properties of TTSs in the Orion Nebula 
Cluster (ONC), which is even younger than the Taurus region. At a distance of 400-500 
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pc ( IHirota et al.ll2007t ISandstrom et al.l 120071 ; IJeffriesI 120071 ). the ONC is one of the most 
extensively studied regions in the sky. Located just in front of the OMC-1 molecular 
cloud, it is the closest site of massive star formation. The age of the ONC is estimated to 
be ~ 1 Myr by com paring the qbsery ed HR diagram with theoretical pre-main sequence 



evolutionary tracks (iHillenbrand 



19971 ). Strong ionizing radiation, primarily from the young 
and bright 06 star 6 1 Ori C, evaporates the surrounding molecular gas, creating a huge 
cavity and exciti ng the famous Orion N ebula. Though much of this region still suffers strong 



extinction (e.g., 



Carpenter et al. 



20011 ) . a few thousand low-mass stars and brown dwarfs 



are visible in the optical. The young stars in the ONC provide an important sample for 
studies such as of the initial mass function of stellar clusters, the properties of circumstellar 
disks, and the evolutionary effects of high-mass stars on low-mass star formation. The 
dynamic star-formation activity in the ONC gives rise to strong stellar jets and outflows, 
result ing in violent int er actions with the ambient molecular environment on various scales 



e.g., 



Bally et al. 



20001 ). 



O'delll ( 1200 ll ) gives a comprehensive review of this region. 



In this paper, we report magnetic-field measurements for 14 TTSs in the ONC that 
allow the first investigation of the magnetic properties of TTSs in this region. In § 2, we 
describe the observations and data reduction procedures. In § 3, we present our analysis 
technique and the magnetic measurements. A discussion of our results is given in § 4, 
followed by our conclusions in § 5. 



OBSERVATIONS 



We conducted our observati ons on UT Decemb er 10-11, 2006, with the Phoenix high- 



resolution near-IR spectrometer (IHinkle et al. 



20031 ) on the 8.1 m Gemini South Telescope 



at Cerro Pachon, Chile (program ID: GS-2006B-C-8). With the K4484 order-sorting filter, 
the grating was oriented to provide wavelength coverage between 2.2190 and 2.2285 jum in 
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air. This configuration allowed the recorded spectra to contain three magnetically sensitive 
Ti I lines as well as less magnetically sensitive lines of Fe I and Sc I. The 0'.'34 (4-pixel) wide 
slit was used to achieve a spectral resolving power of R = X/5X ~ 50,000. For each star, 
two exposures were acquired in nodding mode with an offset of 5" along the 14" long slit. 
Typical seeing was around 0'.'5. Telluric and spectral standards were also observed each 
night. For data reduction, we used custom IDL soft ware developed in Paper I that was 



origin ally written to handle data from the CSHELL (ITokunaga et al 



19901 : 



Greene et al. 



19931 ) spectrometer at the NASA Infrared Telescope Facility. The reduction procedures are 



described in Paper I and required minimal modification to work with the Phoenix data 
format. Seven strong telluric absorption lines were used for wavelength calibration. A 
summary of the ONC star observations is provided in Table [TJ 



3. ANALYSIS AND RESULTS 

We analyze three Ti I lines at 2.2211, 2.2233, and 2.2274 /im using the same procedures 
that are described in detail in Papers I-IV. The effective Lande g-factors of the three lines 
are 2.08, 1.66 and 1.58, respectively. First, for our objects, we obtain estimates of the stellar 
parameters T e g, logg, and wsinz from the literature. Then we construct model atmospheres 
appropriate for the stars and calculate synthetic spectra. We predict the line profiles with 
no magnetic broadening and model the observed excess broadening to measure the surface 
magnetic fields on the TTSs in Orion. In contrast to our analysis of the stars in Papers III 
and IV, we do not have magnetically insensitive CO data to help verify the accuracy of the 
atmospheric parameters. However, two Fe I lines at 2.2257 and 2.2260 fim and one Sc I line 
at 2.2267 /im are available. These three lines have weak sensitivity to magnetic fields with 
effective Lande g-factors of 1.00, 0.74, and 0.50, respectively. So they essentially serve the 
same purpose as the CO lines (i.e., they help ensure that the excess broadening seen in the 
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Table 1: Journal of Observations 



Object 


H.A. (2000) 


Dec. (2000) 


UT Date 


UT Time 


lotal Exposure lime (s) 


2MAbb 05353126-0518559 


05:35:31.26 


-05:24:18.4 


2006 Dec 10 


01:56 


1200 


VT227 On 


05:35:10.58 


-05:21:56.2 


2006 Dec 10 


02:40 


1600 


2MASS 05351281-0520436 


05:35:12.81 


-05:20:43.6 


2006 Dec 10 


03:16 


1800 


V1123 Ori 


05:35:23.44 


-05:10:51.7 


2006 Dec 10 


03:56 


2280 


OV Ori 


05:35:52.76 


-05:12:59.0 


2006 Dec 10 


04:44 


2600 


V1348 On 


05:35:30.42 


-05:25:38.5 


2006 Dec 10 


05:50 


2880 


LO Ori 


05:35:08.22 


-05:37:04.7 


2006 Dec 10 


06:49 


2880 


V568 Ori 


05-35-36 44 


-05-34-1 1 1 


2006 Dec 1 


07:45 


3000 


LW Ori 


05:35:12.20 


-05:30:32.9 


2006 Dec 11 


01:44 


3600 


V1735 Ori 


05:35:22.45 


-05:09:11.0 


2006 Dec 11 


02:53 


3600 


V1568 Ori 


05:35:36.70 


-05:37:41.5 


2006 Dec 11 


04:03 


3600 


2MASS 05361049-0519449 


05:36:10.49 


-05:19:44.9 


2006 Dec 11 


06:35 


3600 


2MASS 05350475-0526380 


05:35:04.75 


-05:26:38.0 


2006 Dec 11 


05:15 


3600 


V1124 Ori 


05:35:47.39 


-05:13:18.5 


2006 Dec 11 


07:46 


3600 



Ti I lines is magnetic in origin). These three lines are included in the fitting. 



In order to construct appropria te model at mospheres, we obtai n stellar atmospheric 



parameters from 



Hillenbrand! (119971 ) . We conv ert 



according to the relation provided in 



Hi 



lenbrandl ( 119971 ) spectral types to T e g 



Johnson! (119661 ) . In a few cases, marked in the T e g 



column in Table [2j we adjust T ef j by 100 or 200 K so that the predicted spectrum better 



mate 
from 



res the Fe I and Sc I lines. To calculate log g, we use the stellar luminosities and masses 



Hillenbrand! (119971 1 and the T e g converted from spectral types. Because of their youth, 



many of our objects have surface gravities lower than 3.5 on t he logarithmic scale, which is 



the lowest available in the NextGen Model Atmospheres grid ( lAUard fe Hauschildt 



1995|). 
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In such cases, in order to avoid extrapolating the model atmosphere grid, we use logg = 
3.5. We expect that doing so produces only a modest (~ 10%) error in our measurements 
of magnetic fields, because our analysis technique is not very sensitive to small errors in T e g 
or logg, as shown by the extensive Monte Carlo analysis in Paper III. The logg values we 
derive from the Hillenbrand paper and those actually used in the analysis are listed along 
with other important stellar quantities in Table [2J 

The magnetic models we use for the ONC stars are constructed in the same way as 
models 1-3 for the TWA stars in Paper IV, where they are fully described. Briefly, the 
stellar surface is treated as being composed of several regions. Each region has a filling 
factor (i.e., a fractional size) and is covered with a single magnetic field strength. The 
sum of all filling factors must be unity. All the models used here include a field-free region 
(5 = 0); however, Papers I and II show that generally the data do not demand the models 
have such a field-free region on the stellar surface. We also assume that these regions 
are spatially well mixed over the surface of the star. Models 1, 2, and 3 allow one, two, 
and three magnetic regions on the stellar surface, respectively, in addition to the field-free 
region. For models 1 and 2, the field strengths are free parameters. For model 3, we have 
field strengths in the three magnetic regions fixed at 2 kG, 4 kG, and 6 kG and fit only 
for the filling factor of each region. The mean magnetic field is just the sum of the field 
strengths (Bi) in the magnetic regions weighted by their corresponding filling factors (/»), 

n 

i.e., B = 2_,Bifi, where n is the number of magnetic regions on the stellar surface in each 
model. 

To compute the model spectrum for a star, we first interpolate the NextGen Model 
Atmosphere grid to the appropriate T e g and log g and construct a model atmosphere specific 
for each star. Then we compute the spectrum for eac h magnetic-field region using the model 



atmosphere and a polarized radiative transfer code (jPiskunov 



19991 ) that assumes a radial 



Table 2: Stellar Parameters 



ID 1 



Object 



Spectral Type T e g 
(K) 



logoff 13 logg c vsini log L M M, 
(kins" 1 ) (L Q ) (M ) 



R, Pro/ AgC° 

(i? ) (days) (10 6 yr) 



867 2MASS 05353126-0518559 CTTS K8 3900 

347 V1227 Ori CTTS K5-K6 4200 

391 2 MASS 05351281-0520436 CTTS Ml 3500* 

731 V1123 Ori NTTS M0/K8 3800 

1021 OVOri CTTS K5-K6 4000* 

850 V1348 Ori CTTS K8 3800* 

321 LO Ori CTTS M0 3900* 

923 V568 Ori CTTS Ml 3664 

381 LW Ori CTTS M0. 5 3800* 

705 V1735 Ori s NTTS K4 4400* 

926 V1568 Ori NTTS K7 4000 

5073 2MASS 05361049-0519449 CTTS K3 4200* 

258 2MASS 05350475-0526380 CTTS M0. 5 3664 

995 V1124 Ori NTTS Ml. 5 3589 



2.98 
3.33 
3.27 
3.20 
3.62 
3.47 
3.49 
3.25 
3.32 
N/A 
3.29 
3.97 
3.34 
3.17 



3.50 
3.50 
3.50 
3.50 
3.62 
3.50 
3.50 
3.50 
3.50 
4.00 
3.50 
4.00 
3.50 
3.50 



11.5±1.0 
10.0±0.8 
10.6±1.5 
16.4±0.8 
14.0±1.8 
15.5±1.3 
12.5±1.7 

6.2±1.9 
15.9±1.2 
16.6±0.8 
15.0±1.0 

9.2±1.2 
10.4±1.0 

6.2±1.3 



0.249 
0.086 
-0.165 
0.098 
0.105 
-0.107 
-0.195 
-0.192 
-0.150 
N/A 
0.111 
0.191 
-0.157 
-0.162 



0.30 
0.41 
0.29 
0.35 
0.63 
0.40 
0.38 
0.29 
0.32 
N/A 
0.40 
1.14 
0.33 
0.26 



2.93 10.66 
1.72 7.33 
2.26 16.33 



2.59 



2.35 



7.64 



2.36 N/A 
2.04 7.78 
1.75 7.74 
2.00 1.15 
1.95 16.20 
N/A N/A 

2.37 5.54 



N/A 



2.08 10.98 
2.15 8.18 



0.7 
3.4 
1.1 
0.8 
1.3 
1.3 
2.0 
1.3 
1.1 
N/A 
1.3 
1.4 
1.2 
1.2 



"The ID, spectral t ype, and Af* are taken from Hillenbrand! (|l997l ). The usini is from 



Sicilia-Aguilar et al 



(|2005|) 



Calculated using Lf, i, M* from 
c Used for this analysis. 
^Rotation periods are taken from 



Hillenbrand! (|1997l ) and T e g used in our models. See text in § 3. 



Flaccomio et al 



(2003) 



Siess et al 



(2000). 



e Stellar ages are calculated from the pre-MS evolutionary tracks of ! 
*T e g is adjusted by 100 or 200 K after converted from spectral type. See text in § 3. 

9 Bolometric luminosity of V1735 Ori is not available in the literature, so there are no estimates of radius, 
surface gravity, mass and age for this star. 

magnetic geometry in the stellar photosphere. We add in rotational broadening to the line 
profiles according to the v sin i of the star and convolve the resulting profiles with a Gaussian 
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function corresponding to the observed spectral resolution (50, 000). The spectrum for each 
region is then weighted by its filling factor and summed to f orm the final model spectrum . 



19921 ) 



We use the nonlinear least-squares technique of Marquardt (IBevington Sz Robinson! 
and solve for the combination of field strengths and filling factors that yields the best fit to 
the observed spectrum. 

The majority of our objects are CTTSs and are subject to substantial K-band veiling. 
Veiling is an excess continuum source, presumably from the circumstellar disk, that acts 
to weaken lines in continuum normalized spectra. Therefore, in addition to solving for 
magnetic field strengths and filling factors, our models also allow K-band veiling (r^-) as a 
free parameter and assume it to be constant across the wavelength range of interest. This 
extra parameter only affects the depths of the lines and does not directly change the line 
widths. The Fe I and Sc I lines help ensure that the veiled Ti I line profiles with little 
magnetic broadening are accurately predicted. The measured mean magnetic field for each 
of the magnetic models (models 1-3), fitted K-band veiling, and corresponding reduced \ 2 
values are listed in Table [3j Typical uncertainty in the veiling measurements is between 
0.05 and 0.10. 

As examples of the fits, the K-band spectra of 2 MASS 05353126-0518559, VI 123 Ori, 
and V1348 Ori are shown in Figures dH31 (The spectra of other targets are presented as 
online-only Figures 7-17.) When correcting the observed spectra for telluric absorption, 
there are only two particularly strong telluric lines that affect the photospheric lines used 
in our analysis. The horizontal bars in Figures HH3] mark the two telluric features that are 
stronger than 3% of the continuum. Though these two regions are li kely to have higher 



uncert ainties than the regions that not affected by telluric absortion, 



Johns-Krull et al. 



( 120091 ) found that arbitrarily increasing the uncertainties of the pixels in Phoenix data that 



are significantly affected by the telluric correction has essentially no effect on the fitting 
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results. The Phoenix data used here is even less affected by telluric contamination and 
we expect uncertainties in its correction to have similarly negligible effect. To show an 
example of how sensitive our analysis is to the uncertainties in T e g, we analyze 2MASS 
05353126-0518559 with a T eff that is ± 300 K from the value we used (3900 K), and the 
measured magnetic field values are changed by only 3% (—300 K) and 8% (+300 K). As 
demonstrated in Paper III, small errors in T e g and log g in general only affect the magnetic 
field measurements by a few percent. 

Additional errors might be incurred by errors in v sinz values, as we do not have 
magnetically insensitive CO lines to check on the atmospheric parameters (in particular 
v sin i) used in the spectrum synthesis. However, the lack of CO line data is not expected to 
be a major issue for this work. In previous papers in this series, the CO data have largely 



served as a check on optically determined v sini measurements. For example, 



Johns- Krull 



(120071 ) did use v sin i measurements determin ed from the CO data 



study. These v sin i values were published in 



Johns-Krull &: Valentil (120011 ) where it was 



or most stars in that 



shown that the CO determined v sin % values are consistent with optically determined values 
such as those used in the present study. To estimate the effect of uncertainties in vsini 
on the recovered values of B, we have repeated the analysis of 2MASS 05353126-0518559 
(medium vsini, strong field) and V568 Ori (low v sini, weaker field) where we have held 
the vsini of each star fixed at values corresponding to ±1 km s" 1 (a typical uncertainty). In 
the case of the 2MASS star, the resulting change in B was 1% and 3%, and in the case of 
V568 Ori the resulting change was 1.5% and 8%. The reason for such a small effect of the 
recovered field is that the magnetic broadening is significantly stronger than the rotational 
broadening in the Ti I lines. To make a final estimate of our field uncertainties from possible 
inaccuracies in stellar parameters, we conservatively assume a 10% uncertainty resulting 
from mismatches in both T e g and logg and add these together in quadrature with an 8% 
uncertainty resulting from errors in the adopted vsini. These sources of uncertainties add 
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up to 16 %. For most of our stars, there is a noticeable decrease in \ 2 f° r models 2 and 3 
relative to the single- component model (1). The two stars for which this is not the case are 
OV Ori and LW Ori, which both have relatively low signal-to-noise spectra. On the other 
hand, models 2 and 3 do not agree perfectly, even though they do generally return mean 
fields that agree to within the 16% uncertainty described above. Therefore, we adopt as 
our final mean field value the average mean field returned by models 2 and 3. We take the 
difference between this mean and the individual values as an estimate of the uncertainty 
of the technique, and add it in quadrature to the 16% uncertainty described above. For 
OV Ori and LW Ori, we average the mean field from all 3 models and take the standard 
deviation as the an estimate of the uncertainty, which is added in quadrature to the 16% 
uncertainty. The final B values and the associated uncertainties are listed in the last two 
columns of Table [3J We use the adopted values of B and its uncertainty, and propagate 
them through all subsequent calculations and comparisons. 



DISCUSSION 



For 14 T Tauri stars in the ONC, we model the Zeeman broadening in three 
magnetically sensitive Ti I lines and measure the mean surface magnetic field strengths 
on these stars. We have attempted to account for various sources of uncertainty in our 
measurements, but there are some potential sources for which this is quite difficult. An 
example of this is the assumed geometry of the field. We have assumed a purely radial 
field geometry in the spectrum synthesis. This was originally motivated by analogy to the 
Sun where observations suggested surface fields (at l east in high flux density regions) are 



primarily oriented along the vertical direction (e.g., 



Stenflo 



20101 ). A number of recent 



Zeeman Doppler Imaging (Z DI) studies have begun t o 



fields on T Tauri stars (e.g. 



Donati et al. 



2007 



2008 



prov i de surface maps of magnetic 



2010 



Hussain et al. 



2009|). These 
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Table 3: Magnetic-Field Measurements. 



Object 


Type 




Model 1 




Model 2 






Model 3 




Adopted 






Xl (B=0) 


B (kG) {xD 


TK 


B (kG) {xD 


TK 


B (kG) (xl) 


TK 


B (kG) 


a (kG) 


2MASS 05353126-0518559 


CTTS 


14.19 


1.91 


(1.99) 


0.85 


2.68 


(1.60) 


1.03 


2.99 


(1.62) 


1.05 


2.84 


0.48 


V1227 Ori 


CTTS 


11.99 


1.81 


(3.00) 


0.72 


2.13 


(2.50) 


0.62 


2.15 


(2.77) 


0.65 


2.14 


0.34 


2MASS 05351281-0520436 


CTTS 


17.85 


1.41 


(1.85) 


0.63 


1.81 


(1.81) 


0.81 


1.58 


(1.82) 


0.66 


1.70 


0.29 


V1123 Ori 


NTTS 


9.71 


1.94 


(3.09) 


0.00 


2.43 


(2.64) 


0.00 


2.58 


(2.69) 


0.00 


2.51 


0.41 


OV Ori 


CTTS 


3.16 


1.55 


(1.10) 


0.59 


1.64 


(1.19) 


0.63 


2.35 


(1.06) 


0.72 


1.85 


0.53 


V1348 Ori 


CTTS 


9.78 


2.05 


(3.54) 


0.38 


3.05 


(3.14) 


0.52 


3.23 


(2.74) 


0.58 


3.14 


0.51 


LO Ori 


CTTS 


9.45 


2.70 


(2.81) 


1.44 


3.26 


(2.52) 


1.61 


3.63 


(2.41) 


1.56 


3.45 


0.58 


V568 Ori 


CTTS 


10.98 


1.32 


(1.83) 


0.64 


1.59 


(1.30) 


0.65 


1.46 


(1.39) 


0.55 


1.53 


0.25 


LW Ori 


CTTS 


8.30 


1.04 


(1.32) 


1.43 


1.32 


(1.33) 


1.56 


1.55 


(1.30) 


1.58 


1.30 


0.33 


V1735 Ori 


NTTS 


3.12 


1.54 


(1.47) 


0.00 


2.14 


(1.38) 


0.00 


2.01 


(1.38) 


0.00 


2.08 


0.34 


V1568 Ori 


NTTS 


4.00 


1.06 


(2.16) 


0.00 


1.42 


(1.77) 


0.00 


1.42 


(2.01) 


0.00 


1.42 


0.23 


2MASS 05361049-0519449 


CTTS 


9.99 


1.06 


(2.10) 


0.24 


2.28 


(1.71) 


0.26 


2.34 


(1.67) 


0.25 


2.31 


0.37 


2MASS 05350475-0526380 


CTTS 


13.95 


1.97 


(2.90) 


0.60 


2.1 SI 


(2.42) 


0.79 


2.76 


(2.44) 


0.70 


2.79 


0.45 


VI 124 Ori 


NTTS 


8.90 


1.45 


(3.90) 


0.00 


2.17 


(3.19) 


0.00 


2.00 


(3.25) 


0.00 


2.09 


0.34 



studies show that the field topology is not purely radial; however, it is not at all clear that 
the ZDI results g ive an accurate picture of the majority of the magnetic flux on the star. 



Reiners fc Basri 



(120091 ) showed that in strongly magnetic M stars, Stokes V based studies 
(ZDI) miss 86% or more of the magnetic flux present compared to Stokes I based studies 
due to the effects of flux cancellation in Stokes V measurements which are not present in 
Stokes I measurements. 

Our work here is based on Stokes I. If the field were entirely aimed along the line 
of sight, or entirely perpendicular to it, the relative strength of the 7r and a Zeeman 
components would change and affect our results somewhat. However, the polarization 



measurements m entioned above show such pathologica 



present (see also 



Johns-Krull et al. 



1999a 



Daouetal 



field orientations are g enerally not 



2006 



Yang et al 



2007|). Instead, 
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2MASS 05353126-0518559 

1.00 

0.95 

| 0.90 

S 0.85 
c 

~Z 2.2210 2.2220 2.2230 

> 

ra 

£ 1.00 
0.95 
0.90 
0.85 

2.2250 2.2260 2.2270 2.2280 

Wavelength(^m) 

Fig. 1. — K-band infrared spectra of 2MASS 05353126-0518559 are shown in the histogram. 
The spectral regions shown in bold are the magnetically sensitive Ti I lines actually used in 
the fit. The dash-dotted lines show a model with no magnetic field. The smooth lines are 
the best fit with magnetic broadening based on model 3. The solid horizontal bars mark the 
wavelengths where telluric absorption is stronger than 3% of the continuum and coincides 
with the lines used in the fit. 





there is likely a stong mixture of field orientations on the star: some parallel to the line of 
sight, some perpendicular to the line of sight, and all angles in between. While we assume a 
radial geometry on the star, when we integrate over the visible stellar disk, we then include 
a good mixture of field lines with all orientations to the line of sight. As a result, we do 
not expect our results to be strongly affected by our assumption of a radial field geometry 
Again, as mentioned above, sin ce the ZDI results like ly do not provide good guidance for 



the geometry of the total field (IReiners fc Basri 



20091 ). this is probably the best we can do 



until better guidance comes from theory or even better observations. 

For most of our targets, the single magnetic component fit of model 1 is not good 
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V1123 Ori 

1.05 
1.00 
0.95 



1 0^80 

£ 0.75 
c 

o 2.2210 2.2220 2.2230 

> 

« 1.05 
-2 1.00 

0.95 

0.90 

0.85 

0.80 

0.75 

2.2250 2.2260 2.2270 2.2280 

Wavelength(^im) 



Fig. 2. — K-band infrared spectra of V1123 Ori are shown in the histogram. The spectral 
regions shown in bold are the magnetically sensitive Ti I lines actually used in the fit. The 
magnetically sensitive Ti I lines are in bold. The dash-dotted lines show a model with no 
magnetic field. The smooth lines are the best fit with magnetic broadening based on model 
3. The solid horizontal bars mark the wavelengths where telluric absorption is stronger than 
3% of the continuum and coincides with the lines used in the fit. 





enough to account for all the excess broadening in the observed line profiles. For a few 
targets such as OV Ori and LW Ori, adding more magnetic components does not improve 
the fits significantly. That is due to a combination of low signal-to-noise ratios and small 
magnetic broadening. Models 2 and 3 generally yield higher mean magnetic field strengths 
and better matches to the observations, resulting in lower reduced x 2 values. In most cases, 
the mean field values recovered by models 2 and 3 are in good agreement with similar 
reduced x 2 values. On average, the difference between the two recovered values is less than 
9%. 

Since we apply models that are consistent with previous studies of 14 stars in Taurus 
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V1 348 Ori 



c 






> 




2.2210 



2.2220 



2.2230 




2.2250 2.2260 2.2270 

Wavelength(i^m) 



2.2280 



Fig. 3. — K-band infrared spectra of V1348 Ori are shown in the histogram. The spectral 
regions shown in bold are the magnetically sensitive Ti I lines actually used in the fit. The 
dash-dotted lines show a model with no magnetic field. The smooth lines are the best fit 
with magnetic broadening based on model 3. The solid horizontal bars mark the wavelengths 
where telluric absorption is stronger than 3% of the continuum and coincides with the lines 
used in the fit. The emission feature at 2.2227 /zm is a molecular hydrogen line. 



( 1 Johns- Krull 



20071 ) and 5 stars in the TW Hydrae association (Papers III and IV), we are 
able to directly compare the measurements of all 33 stars in the three regions. Out of the 33 
stars in total, 27 stars have spectral types between K7 and M3. The rest of the sample have 
spectral types between K3 and K6, except for the K0 star T Tau. Therefore, these stars 
possess similar spectral properties and allow a meaningfu l comp arison. We calculate stellar 



ages using the pre-MS evolutionary tracks of ISiess et al. 



( 120001 ). and plot the measured 



magnetic field values and magnetic fluxes (&irR* 2 B) against stellar ages, shown in Figures 
H] and [5j respectively. We note that while stellar magnetic field strength spans a range of 
values and displays no apparent correlation with stellar age, stellar magnetic flux decreases 
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consistently with time. The primary reason for the inverse correlation of magnetic flux with 
age is the diminishing stellar radii as the stars evolve down their Hyashi tracks. However, as 
shown below, the measured fields/fluxes do not correlate with any other expected dynamo 
properties, including for example the stellar gravity. As a result, we find no evidence the 
the field strengths are being set by some other property of the stars, which might produce 
an inverse correlation with age/radius as a by-product. The observed trend of the flux 
decreasing with age supports the i dea that the strong magnetic fields of TTSs could be of 



primordial origin. As suggested by 



Mossl (120031 ) . during the early stages of stellar evolution, 



while primordial magnetic flux does decay, young stars may be able to maintain part of the 
fossil flux from the molecular cloud material, especially later when a raditive core is formed. 
It should be noticed that currently there are only a few measurements for stars younger 
than 1 Myr or older than 10 Myr, so more observations are needed to further trace out the 
evolution of stellar magnetic fields with time. 

In addition to a fossil origin for the magnetic fields on young stars, these observed 



strong fie 



ds could in principle be the result o f dynamo action, either a solar-like a — Q 



dynamo (jDurney &: Latour 



convective dynamo (e.g., 



1978 



Parker 



Durney et al. 



19931) or as a result of a distributed, turbulent or 



19931 ) . In current a — Q dynamo models, magnetic 



fields are generated by differential rotation at the region between the radiative interior 
and the outer convective shell, a region known as the tachocline. In such a dynamo, it 
is believed that levels of magnetic activity are related to the Rossby number, defined as 
Prnt/Tr., whe r e P rn t is the rotation period and r c is the convective turnover time. However, 



Johns-Krulll (120071 ) found no correlation for 14 Taurus stars and TW Hya between the 
measured magnetic field values a nd various stellar o r dynamo paramete rs. Here, for a 



combined sample of 31 stars from 



Yang et al. 



(|2008|) 



Johns-Krulll (120071 1 and this study. 



we look for a correlation between measured magnetic prope rties and various potenti al 



dynamo parameters. We obtain r c from models described in 



Kim &: Demarqud (119961 ) and 



1.0 



1 



10 

Age (Myr) 



Fig. 4. — Measured magnetic fields plotted against age for the Orion, Taurus/Auriga and 
TWA Stars. The triangles represent the Orion stars from this work. The filled circles 
represent the TWA s tars from Papers H I and IV, and the hollow diamonds are the Tau- 



rus/Auriga stars from I Johns-Krulll (120071 ). The sizes of error bars are 10% for the TWA and 



Taurus/Auriga stars and the values in the last column of Table [3] for the Orion stars. 



Preibisch et al. 



( 120051 ). In Table HI we give the Pearson linear correlation coefficient, r, 
and its false alarm probability, f p , for several comparisons. No significant correlation is 
seen between the observed magnetic field strength/flux and any dynamo parameter. Such 
lack of correlation is generally contrary to the prediction of theories of solar-like interface 
dyna mos. Besides interfa ce dynamos, there are models that involve a distributed dynamo 



(e.g. 



Durney et al. 



19931 ). This kind of dynamo is thought to be able to operate in fully 



convective stars such as TTSs and create large-scale magnetic fields via convective motions. 
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Fig. 5. — Measured magnetic flux plotted against age for the Orion, Taurus/Auriga and 
TWA Stars. The triangles represent the Orion stars from this work. The filled circles repre- 
sent the T WA stars from Pap ers III and IV, and the hollow diamonds are the Taurus/Auriga 



stars from 



Johns-Krulll (l2007t ). 



Recently, 



Christensen et al 



( 120091 ) extended a scaling law derived from convection-driven 
geodynamo models to the regime of stars, including TTSs, and their predictions for the field 
strengths on TTSs are generally consistent with the observed values. However, magnetic 
fields generated by a distri buted dynamo are also expected to haye a d ependence, though 



possibly weak, on rotationt lDurney et al.lll993t IChabrier fc Kiiker 



20061 ). which we do not 



observe. Therefore, we suggest that the survival of the fossil magnetic flux seems to be a 
more promising explanation for the origin of strong magnetic fields on TTSs; however, more 
work in this area is needed. 
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Table 4: Correlations of Measured Magnetic Parameters with Dynamo Parameters 



Quantities Compared 


r 


fp 


B vs. P'l 


-0.01 


0.94 


$ vs. p- a \ 


-0.04 


0.83 


B VS. T c 


0.24 


0.20 


$ VS. T c 


-0.08 


0.64 


B vs. t c P~1 


0.05 


0.79 


$ VS. T c P r 7 a \ 


-0.15 


0.42 


B vs. log g 


0.14 


0.46 


$ vs. log g 


-0.45 


0.01 



We also examine the relationshi p between magnetic flux and the X-ray properties 



of the ONC stars. As first noted by 



Pevtsov et al. 



(120031 ) for the Sun and a number of 



main-sequence dwarfs, there is a nearly linear correlation between X-ray luminosities, L x , 
and the total unsigned magnetic flux, $. The correlation holds over 12 orders of magnitude 
in both quantities. We use this relationship with our measured magnetic flux values to 



calculate the predicted X-ray luminosities for our ONC stars. These va 
able [51 w here we compa r e them with the observed values reported by 



ues are listed in 



Flaccomio et al. 



(120031) and 



Getman et al. 



(120051 ). The typical uncertainty in the observed L x is 0.2 dex. 
The comparison is also plotted in Figure El On average, the predicted L x values of the 
ONC stars are larger than the observed ones by a factor of 4.9 ± 1.4. This suppression of 
the X-ray emission is consistently found as well for the Taurus and TWA stars, also shown 
in Figure EJ though the reduction is different in each region. 
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Large flares are thought to be responsible for harder X-ray emission on TTSs, while 
softer, persistent X-ray emission may be generated by so-called mic roflares that relea se 



energy from stress ed magnetic fields in the coronae of these stars ( jGiidel et al. 



Arzner et al. 



20031 : 



20071 ). On the Sun, these stresses are believed to be built up by convective 
motions acting on the footpoints of coronal loops. For TTSs, strong surface magnetic fields 
may decrease the efficiency by which convective gas motions are able to move footpoints 
around and build up these magnetic stresses. Thus, in addition to listing the measured 
magnetic field strengths based on model 3, B, we also list in Table |5] the equipartition field 
strength, B eq = (87rP g )^ 2 . This quantity is calculated assuming that the magnetic pressure, 
B eq 2 /(8n), is in balance with the surrounding unmagnetized photospheric gas pressure, P g . 
We use the gas pressure from a location in the model atmospheres where the local gas 
temperature is approximately equal to the effective temperature of the star. This level in 
the stellar photosphere is in the vinicity of where the observed continuum forms, so the gas 
pressure here should be generally greater than the regions where the Ti I lines form, setting 
an upper limit. The ratio of the observed field strength to the equipartition field strength 
ranges from 1.3 to 3.9 for the ONC stars, indicating the dominance of magnetic pressure 
over gas pressure in the photospheres of these TTSs. This dominance of magnetic pressure 
may partially prohibit the production of X-ray emission. Since these ratios of the ONC 
stars are close to those of the TWA stars, it is likely no coincidence that in Figure [6] the 
observed X-ray luminosities of most ONC stars are smaller than the predicted values by a 
similar factor as those of the TWA stars. 



5. Conclusion 

By modeling the Zeeman broadening in three magnetically sensitive Ti I lines, we 
consistently measure kilogauss-level magnetic fields in the photospheres of 14 TTSs in 
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Table 5: Magnetic-Field Properties and X-Ray Luminosities. 



Object 


B 


eq 


Predicted log Lx 


Observed logL^ 




(kG) 


(kG) 


(ere/s) 


(ere/s) 


2MASS 05353126-0518559 


2.84 


0.99 


31.16 


30.37 


V1227 Ori 


2.14 


0.87 


30.49 


30.20 


2MASS 05351281-0520436 


1.70 


1.03 


30.65 


29.65 


V1123 Ori 


2.51 


1.00 


30.98 


30.65 


OV Ori 


1.85 


1.07 


30.73 


N/A 


V1348 Ori 


3.14 


1.00 


30.86 


30.63 


LO Ori 


3.45 


0.99 


30.75 


<30.23 


V568 Ori 


1.53 


1.07 


30.47 


30.04 


LW Ori 


1.30 


1.00 


30.37 


29.63 


V1735 Ori 


2.08 


1.35 


N/A 


N/A 


V1568 Ori 


1.42 


0.93 


30.61 


<30.38 


2MASS 05361049-0519449 


2.31 


1.42 


30.85 


30.69 


2MASS 05350475-0526380 


2.79 


1.02 


30.81 


30.45 


V1124 Ori 


2.09 


1.67 


30.70 


30.29 



the ONC. We find a systematic decrease of stellar magnetic flux from the young Orion 
region (~ 1 Myr), through the Taurus region (~ 2 Myr), to the TWA (~ 10 Mys). 
This finding supports a primordial origin of the magnetic fields on TTSs, though recent 



convective dynamos by 



Christensen et al. 



(120091 ) are able to produce magnetic fields that 



have strengths consistent with observations. Because convective motions are likely partially 
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Fig. 6. — Observed X-ray luminosity plotted against the predicted X-ray luminosity for the 
Orion, Taurus/Auriga and TWA stars. The solid line is the line of equality, and the dashed 
line is where the predicted value is 10 times the observed one. The triangles represent the 
Orion stars from this work. The filled circles represent the TWA stars fro m Paper IV, and 



the hollow diamonds are the Taurus/Auriga stars from iJohns-Kruill (120071 ). 



inhibited by the dominance of magnetic pressure over unmagnetized gas pressure in the 
stellar photospheres, the strong magnetic fields could also be responsible for suppression of 
X-ray emission on TTSs relative to that expected from main sequence star calibrations. 



This paper is based on observations obtained with the Phoenix infrared spectrograph, 
developed and operated by the National Optical Astronomy Observatory. 
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A. Online-Only Figures 

V1127 Ori 
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Fig. 7. — K-band infrared spectra of V1227 Ori are shown in the histogram. The spectral 
regions shown in bold are the magnetically sensitive Ti I lines actually used in the fit. The 
dash-dotted lines show a model with no magnetic field. The smooth lines are the best fit 
with magnetic broadening based on model 3. The solid horizontal bars mark the wavelengths 
where telluric absorption is stronger than 3% of the continuum and coincides with the lines 
used in the fit. 
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2MASS 05351 281 -0520436 
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Fig. 8.— K-band infrared spectra of 2MASS 05351281-0520436. 
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Fig. 9. — K-band infrared spectra of OV Ori. 
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Fig. 10. — K-band infrared spectra of LO Ori. 




Fig. 11. — K-band infrared spectra of V568 Ori. 
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Fig. 12. — K-band infrared spectra of LW Ori. 




Fig. 13.- 



K-band infrared spectra of VI 735 Ori. 
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Fig. 14. — K-band infrared spectra of V1568 Ori. 




Fig. 15.- 



K-band infrared spectra of 2MASS 05361049-0519449. 
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Fig. 16.— K-band infrared spectra of 2MASS 05350475-0526380. 




Fig. 17.- 



K-band infrared spectra of VI 124 Ori. 
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